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Abstract 

Large-distance modification of gravity may be the meclianism for solving 
the cosmological constant problem. A simple model of the large-distance 
modification — four-dimensional (4D) gravity with the hard mass term — 
is problematic from the theoretical standpoint. Here we discuss a different 
model, the brane-induced gravity, that effectively introduces a soft graviton 
mass. We study the issues of unitarity, analyticity and causality in this model 
in more than five dimensions. We show that a consistent prescription for 
the poles of the Green's function can be specified so that 4D unitarity is 
preserved. However, in certain instances 4D analyticity cannot be maintained 
when theory becomes higher dimensional. As a result, one has to sacrifice 
4D causality at distances of the order of the present-day Hubble scale. This 
is a welcome feature for solving the cosmological constant problem, as was 
recently argued in the literature. We also show that, unlike the 4D massive 
gravity, the model has no strong-coupling problem at intermediate scales. 
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1 Large distance modification of gravity: 
formulating the problem 

The reason underlying the observed acceleration of the universe is puzzling. It 
could be a tiny amount of vacuum energy. However, this possibility is hard to 
reconcile with known particle-physics models. Instead, it might well be that a new 
physical scale exists in the gravitational sector and the laws of gravity and cosmology 
are modified at this scale. To be consistent with data and be able to predict the 
accelerated expansion, the new scale should be roughly equal to Hq^ ~ 10^^ cm — 
the present-day value of the Hubble length. In this regard, developing models in 
which gravity gets modified at cosmological distances, becomes a timely endeavor. 
A generally covariant theory of the large-distance modification of gravity is the DGP 
model [1]. The gravity action of the model can be written as follows: 

S J d'x^R{g) +^— y d'xd^'y^gn.Mg), (1) 

where R and TZ^+n are the four- dimensional and (4 -|- A'^)-dimensional Ricci scalars, 
respectively, and stands for the gravitational scale of the bulk theory. Extra 
dimensions are not compactified, they asymptote at infinity to Minkowski space. 
The higher-dimensional and four-dimensional metric tensors are related as 

g{x,y = Q) = g{x). (2) 

The first term on the right hand side of (1) acts as a kinetic term for a 4D graviton 
while the second term acts as a gauge invariant mass term. The observable matter 
is assumed to be localized on a 4D surface y = 0. 

The present work is devoted to the study of the DGP scenario in the case N > 2 
(see Ref. [2], [3]). Such models have a string theory realization [4]. More impor- 
tantly, these models are potential candidates for solving [5, 6] the cosmological con- 
stant problem (see also Refs. [7]-[24] for interesting cosmological and astrophysical 
studies). 

The equation of motion for the theory described by the action (1) takes the form 

5W(y) 5^ 5^ + M^^ G(g = -T,^{x) S'^ ^^(y) . (3) 

Our conventions are as follows: 

TjAB = diag [+--...-], A,B = 0,1,...,3 + N, 

= 0,1,2,3, a,6 = 4,5, ...,3 + A^. (4) 

G^^}, and denote the four-dimensional and D-dimensional Einstein tensors, re- 
spectively. We choose (for simphcity) a source locahzed on the brane, T^jj{x)5^^\y). 
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Gravitational dynamics encoded in Eq. (3) can be inferred both from the four- 
dimensional (4D) as well as (4 + iV) -dimensional standpoints. From the 4D perspec- 
tive, gravity on the brane is mediated by an infinite number of the Kaluza-Klein 
(KK) modes that have no mass gap. Under conventional circumstances (i.e., with 
no brane kinetic term) this would lead to higher- dimensional interactions. However, 
the large 4D Einstein-Hilbert (EH) term suppresses the wave functions of heavier 
KK modes, so that in effect they do not participate in the gravitational interac- 
tions on the brane at observable distances [25] . Only light KK modes, with masses 
rriKK ^ rUc, 

remain essential, and they collectively act as an effective 4D graviton with a typical 
mass of the order of rric and a certain smaller width. 

Assuming that ~ 10~^ eV or so, we obtain rUc ^ Hq 10~^^ GeV. Therefore, 
the DGP model with N >2 predicts [26] a modification of gravity at short distances 
M^^ ~ 0.1 mm and at large distances ~ Hq^ ~ 10^^ cm, give or take an or- 
der of magnitude. Since gravitational interactions, nevertheless, are mediated by 
an infinite number of states at arbitrarily low energy scale, the effective theory (1) 
presents, from the 4D standpoint, a non-/oca/ theory [5]. Moreover, as was suggested 
in [27], nonlocalities postulated in pure 4D theory can solve an "old" cosmological 
constant problem [27], and give rise to new mechanisms for the present-day accel- 
eration of the universe [27, 28]. (It is interesting to note that the nonlocalities in 
a gravitational theory that are needed to solve the cosmological constant problem 
could appear from quantum gravity [29] or matter loops [30] in a purely 4D context.) 

On the other hand, from the (4 + iV) -dimensional perspective, gravitational 
interactions are mediated by a single higher-dimensional graviton. This graviton 
has two kinetic terms given in Eq. (1), and, therefore, can propagate differently on 
and off the brane. Namely, at short distances, i.e. at r < ~ H^^ ~ 10^^ cm, 
the graviton emitted along the brane essentially propagates along the brane and 
mediates 4D interactions. However, at larger distances, the extra-dimensional effects 
take over and gravity becomes (4 + A^)-dimensional. 

As was first argued in Ref. [2], the results in > 2 DGP models are sensitive 
to ultraviolet (UV) physics, in contradistinction to the A^ = 1 model [1]. In other 
words, one should either consistently smooth out the width of the brane [26], or 
introduce a manifest UV cutoff in the theory [31, 32, 26], or do both. With a 
finite thickness, more localized operators appear on the worldvolume of the brane, 
in addition to the worldvolume Einstein-Hilbert term already present in Eq. (1) [26]. 
For instance, one could think of a higher- dimensional Ricci scalar smoothly spread 
over the worldvolume [3]. 

In general, terms that are square of the extrinsic curvature can also emerge. 
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Some of these terms can survive in the hmit when the brane thickness tends to zero 
(i.e. in the low-energy approximation). For instance, in the zero-thickness hmit of 
the brane the following terms might be important: 

5^''\y)K^lK. 5^''\y)h^''d^d,K, 5^''\y)h^^dad,h^\ (6) 

where h denotes small perturbations on fiat space. Although the main features 
of the model, such as interpolation between the 4D power-law behavior of a non- 
relativistic potential at short distances and the higher-dimensional behavior at large 
distances, are not expected to be changed by adding these terms, nevertheless, the 
tensorial structure of a propagator could in general depend on these terms and self- 
consistency of the theory may require some of these terms to be present in the 
actions in a reparametrization invariant way. 

In the low-energy approximation the exact form of these "extra" terms and 
their coefficients are ambiguous, because of their UV origin. They will be fixed in 
a fundamental theory from which the DGP model can be derived [4, 33]. In the 
present paper, in the absence of such a fundamental theory, (but in the anticipation 
of its advent), we would like to study a particular parametrization of these "extra" 
terms, for demonstrational purposes. According to our expectations, physics in the 
self-consistent theory will have properties very similar to those discussed below. 
We will show that these properties are rather attractive since they do avoid severe 
problems of 4D massive gravity. 

Consider the action 

S ^ ^ J d'x^ (a Rig) + bU^+N) + — ^ / d'x d^'y V^TZ^+Nig) , (7) 

where, in addition to the 4D EH term, a D-dimensional EH term localized on the 
brane is present. Here a and b are some numerical coefficients. We will study 
the properties of the system described by (7) for different values of a and b. The 
action (7) is fully consistent with the philosophy of Ref. [1]: if there is a (1+3)- 
dimensional brane in D-dimensional space, with some "matter fields" confined to 
this brane, quantum loops of the confined matter will induce all possible structures 
consistent with the geometry of the problem, i.e. (l-|-3)-dimensional wall embedded 
in D-dimensional space. 

The equation of motion in the model (7) takes the form 

5W(y)M|i (aGg] + 51(5^ + M^^gS = -T,^{x) S>^d^^d^''\y) . (8) 

In deriving the above equation we first introduced a finite brane width A, and then 
took the A — * limit in such a way that no surface terms appear. In general, the 
results depend on the regularization procedure for the brane width. In the present 
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work we adopt a simple prescription in which derivatives with respect to the trans- 
verse coordinates calculated on the brane vanish in the A — > limit (a unique 
prescription could only be specified by a fundamental theory.) . As previously, G^^^ 
and G^^^ denote the four-dimensional and D-dimensional Einstein tensors, respec- 
tively, while a and b are certain constants. In order to be able to describe 4D gravity 
at short distances with the right value of the Newton's coupling we set 

a + 6 = 1 . (9) 

Note, that the first two terms in parenthesis on the left-hand side of Eq. (8) can be 
identically rewritten as 

(a + 6)G(fj +b{-d,d,K - d.dah^ - + r),,dlh% 

+ d^d.K - Vf.udl K + 2r/^.a„a«/i'^" + v^Mh'^'') . (lo) 

The above equation of motion (8) - which should be viewed as a regularized version 
of the DGP model - could be obtained from the action (1) as well, provided the 
latter is amended by certain extrinsic curvature terms (for more details see Ref. [33]). 
Below we will study this version of the regularized DGP model for various values of 
the parameters a and h. Certain issues in the a — h — 1 case, regularized with 
finite brane width, have been recently analyzed in Ref. [3]. We will find below that 
phenomenologically more attractive is the a — 1, h — Q case. 
The issues to be addressed are: 

Assume that gravity measurements are done at points xi and X2 that are confined 
to the brane. At distances |xi — X2I < Hq^ the graviton propagator imitates that 
of a massive 4D unstable particle with mass ~ rric- Given that the model itself is 
intrinsically (4-|-A'")-dimensional, the following questions must be answered: 

(i) Does the graviton propagator ^2, y = 0) satisfy the requirements of 
four-dimensional unitarity? 

(ii) Do "abnormalities" occur at 4D momenta much smaller than the ultraviolet 
(UV) cut-off and larger than the infrared (IR) crossover scale, such as a precocious 
onset of the strong-coupling regime? 

Needless to say, the answers to these questions determine whether the DGP 
model is intrinsically self-consistent and phenomenologically viable. The answer to 
the first question will be demonstrated to be positive while to the second negative. 
That is to say, the situation is most favorable. We hasten to add that it is not 
trivial to see that this is indeed the case. It is necessary to carry out a rather subtle 
analysis which circumvents stereotypes in, at least, one point. 

For what follows it is instructive to confront the DGP model with the 4D "hard"- 
massive gravity [34] (or the Pauli-Fierz (PF) gravity) which leads to a large-distance 
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modification of interactions too. (Tlie action of the "hard" -massive gravity is given 
below in Eq. (98).) In particular, we compare the perturbative treatments of these 
two models. Perturbation theory in Newton's constant in the 4D massive gravity 
breaks down at a scale much lower than the cut-off scale of the theory — this was first 
obtained for spherically symmetric sources in Ref. [35]. The origin of this breaking 
can be traced back [36] to Fcynman diagrams involving nonlinear interactions of 
gravitons. In terms of degrees of freedom, it is the longitudinal polarizations of the 
massive gravitons that are responsible for the perturbation theory breaking. This 
can be readily inferred from dynamics of these modes analyzed in Ref. [37] . (Note 
also that the PF graviton propagates six degrees of freedom instead of five [38] in 
the full quantum theory. This makes the corresponding Hamiltonian unbound from 
below [38]. As a result, solutions exist that destabilize the empty Minkowski space 
[39] and the instability can develop practically instantaneously.). 

As will be shown, in this respect the DGP gravity (7) presents a drastic im- 
provement. In contradistinction with the 4D Pauli-Fierz theory of massive gravity, 
the precocious breakdown of perturbation theory does not occur in the DGP model 
with N > 2. A direct analogy with the Higgs mechanism for non-Abelian gauge 
fields is in order here. 

For non-Abehan gauge fields with the hard gauge-boson mass term, appropriate 
nonlinear amplitudes invalidate the perturbative expansion (i.e. violate the unitar- 
ity bound) at a scale set by the gauge-boson mass divided by the gauge coupling 
constant. To cure this disaster, one introduces an extra scalar — the Higgs field. ^ 
By the same token, certain nonlinear perturbative amplitudes of the 4D "hard"- 
massive gravity blow up precociously [35, 36, 37] at a scale significantly lower than 
a naive UV cut-off of the theory under consideration. The unwanted growth of the 
amplitudes is canceled, however, in the DGP model at the expense of introducing 
an infinite number of the KK fields. Thus, the action (7) gives rise to a gravitational 
analog of the Higgs mechanism, with an infinite number of the "Higgs" fields which 
include both vector and scalar states. 

It has been recently argued [40] that the spectrum of the DGP model contains 
tachyonic states with a negative norm ("tachyonic ghosts"). The conclusion was 
based on an analysis of the poles in the graviton propagator derived from the action 
(1). In fact, the analysis of Ref. [40] leaves aside subtle points of appropriately 
defining the Green's function poles. We formulate and discuss an appropriate rule 
for defining the would-be poles. With this rule accepted, the 4D Tinitarity of the 
Green's functions is guaranteed. This is not the case with respect to 4D analyticity 
and causality, though. In certain instances we will have to sacrifice 4D causality at 
distances of the order of today's Hubble scale. As was argued in Ref. [27] this is a 
welcome feature for solving the cosmological constant problem. 

^The mass of the Higgs field itself needs a stabilization mechanism. This is a separate story, 
however. 
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The organization of the paper is as follows. In Sect. 2 we discuss in detail 
a simplified version of the phenomenon, with the tensorial structure suppressed. 
We consider a scalar field with the Lagrangian similar to Eq. (1) and derive the 
Green's function. There are complex poles of the Green's function on the second and 
subsequent non-physical Ricmann sheets. This corresponds to a resonance nature 
of the 4D massive scalar. In Sect. 3 we discuss the same problem for gravity, i.e. 
including the graviton tensorial structure. The issue of proper definition of poles 
in the graviton propagator emerges in earnest in the trace part. In certain cases, a 
naive way of defining the poles leads to inconsistent results — violations of unitarity. 
For these cases we point out a way out, treating the would-be poles in the Green's 
function in such a way that unitarity is not violated. Finally, in Sect. 4 the N > 2 
DGP model is argued to have no strong coupling problem at intermediate scales, in 
contradistinction to the 4D Pauli-Fierz gravity with the hard mass term. 

2 The simplest example: scalar field 

To warm up, we start our discussion with a simple model of a scalar field $ in (4-l-A'")- 
dimensional space-time. For convenience we separate the dependence of the scalar 
field $ on four- dimensional and higher-dimensional coordinates, Ua) = ^{x, y). 

The two kinetic term action — the scalar counterpart of Eq. (1) — has the form 

S =^ j ci'a;a^$(x,0)9'^$(a:,0) + — ^ / d^xd''ydA^{x,y)d^^{x,y). (11) 

It is important to understand that in the scalar case the analog of the new term 
included in Eq. (8) but absent in (3) reduces, identically, to the already existing 
localized term. This is a consequence of our choice of the regularization of the brane 
width A and the boundary conditions according to which transverse derivatives 
vanish on the brane in the A — > limit. 

To study interactions mediated by the scalar field we assume that $ couples to 
a source J localized in the 4D subspace in a conventional way, / 0) J{x). 

Then the equation of motion takes the form 

(5W(y)M|ia2$(x,0) + Mr^ai$(a;,y) = J(a;)(5W(y). (12) 

The very same equation applies to the scalar field Green's function. 

2.1 Spectral representation 

At first, using the scalar field example, we will summarize general arguments for 
the existence/ absence of a spectral representation in higher-dimension theories with 
the worldvolume kinetic terms as in Eq. (11). Explicit formulae below refer to 
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the scalar case. In the next section we will consider gravity, with the appropriate 
tensorial structure, and will emphasize crucial differences between the present scalar 
example and full-blown gravity. 

By the spectral representation we mean the Kallen-Lehmann (KL) representa- 
tion for the free tree-level propagator in the model (11) in terms of four- dimensional 
Mandelstam variables. Since the theory described by (11) is intrinsically higher- 
dimensional, it is not clear a priori why the spectral representation in terms of the 
4D variables should hold at all. Indeed, on the one hand, the KL representation 
expresses the fact that a given amplitude, as a function of (the 4D momentum 
squared), is analytic in the complex plane everywhere except possible isolated 
poles plus a branch cut along the real positive semiaxis. On the other hand, the 
model (11) is nonlocal from the 4D standpoint; hence, it is not obvious why analyt- 
icity of the amphtude with respect to the 4D variable should take place in the 
regime where higher-dimensional effects become crucial. 

There is an alternative point of view on the KL representation of the model (11). 
We can assume that the extra dimensions y are compactified, with a finite (albeit 
arbitrarily large) compactification radius R. Then, the spectrum of the theory must 
consist of "discretized" Kaluza-Klein modes. Prom the 4D standpoint they are 
just certain massive states. Then one could certainly obtain the KL representation 
by writing the (tree-level) 4D propagator summing up the entire tower of the KK 
eigenstates. 

This strategy is readily implemented in the conventional compactifications, when 
the brane worldvolume term is absent. In this case the spectrum of the Kaluza-Klein 
eigenmodes, 

^x,y)=Y,^n{x)My), (13) 

n 

becomes trivially discrete, with the eigenvalues ~ v? /B?. If so, the expression 
for the Green's function G{p, y — 0) takes the form 

G(p,, = o) = i: ■ p' = P.r^- (14) 

n mi - p^ - le 

In other words, G{p, y = 0) is the sum over the infinite number of poles, with the 
positive-definite residues. As i? — > oo the sum goes into the standard dispersion 
integral, 

G{p,y^O)^- r dt . , (15) 

TT JO t — p^ — le 

where p{t) is a positive-definite spectral density. 

The argument above, as well as the simple representation (14) or (15) following 
from it, neglects the existence of the brane worldvolume kinetic term in the action 
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(the first term on the right-hand side of (11)). This term is crucial, and by no means 
can be neglected. It gives rise to kinetic mixings of the KK modes on the brane 
worldvolume, 

Therefore, the KK modes defined in Eq. (13) are not the eigenstates of the Hamil- 
tonian in the presence of the brane kinetic term. Diagonalization is needed. For 
the scalar field example (11), explicit diagonalization is possible and was in fact 
carried out [25]. As a result, the spectral representation can be argued to exist in 
the desired form, Eq. (15). 

In the case of gravity things are more complicated, however. The worldvolume 
EH term gives rise to kinetic mixings between the massive KK modes of distinct 
spins. It is not obvious how to diagonalize the full linearized Hamiltonian. Even if 
the diagonalization is possible it is not clear whether the diagonal eigenstates are 
states of a definite 4D spin, and not the mixed states. This is all because of the 
large kinetic mixings between all the KK states on the brane. Thus, the spectral 
representation we look for is hard (if possible at all) to obtain through explicit 
summation of the eigenstates of the Hamiltonian. The best one can do is to write 
down the spectral representations in the limiting regimes when the 4D EH term is 
either dominant or negligibly small. We will return to this issue in the Sect. 3. Prior 
to delving in the gravity problem we want to conduct detailed studies of the scalar 
example (11). 

2.2 Solving Eq. (12) in the general case 

To solve this equation it is convenient to Fourier-transform it with respect to "our" 
four space-time coordinates p^, keeping the extra y coordinates intact. Mark- 

ing the Fourier-transformed quantities by tilde, 

$(x,y) ^ l>(p,y), (17) 

we then get from Eq. (12) 

5(^)(t/) M|,(-/) ^(P, 0) + M^^ (V - A,)l>(p, y) = J{p) (5(^)(y) , (18) 
where = Pq — Pi — P2 ~ Ph notation 

- E 5-i (19) 

a=l ^Ua 

is used. 
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We will look for the solution of Eq. (18) in the following form: 

^p,y) = D{p,y)x{p): (20) 

where the function D is defined as a solution of the equation 

{-p^-Ay-ie)D{p,y) = ^^(y) . (21) 

Note that the function D is uniquely determined only after the ie prescription spec- 
ified above is implemented. We also introduce a convenient abbreviation 

Do{p) = D{p,y = 0). (22) 

Now, it is quite obvious that a formal solution of Eq. (18) can be written in 
terms of the function D as follows: 

where $hom(p, y) is a general solution of the corresponding homogeneous equation 
(i.e., Eq. (18) with the vanishing right-hand-side), and c is an arbitrary constant. 
Equation (23) presents, in fact, the Green's function too, up to the factor J(p)/Mpj, 
which must be amputated. In particular, for the Green's function on the brane we 
have 

Gip, 0) - ^ ^p, y = 0) , Ghom(p, 0) = ^ $hom(p, y = 0) , (24) 
J (p) J [p) 

while for arbitrary values of y 

Gip,y) = - 2nl^\^^ N + cG^UP,y), (25) 
p^ Dq(p) - u'"' 



where 



N 



The presence/absence of the homogeneous part is regulated by the ie prescription. 
Note that if the first term on the right-hand side of Eq. (25) has poles on the real 
axis of p^, then the homogeneous equation has a solution 

Ghom(p, y) = D{p, y) 5 (p'Doip) - u"") . (27) 

This fact will play an essential role for gravity, as will be discussed in due course in 
Sect. 3. 
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In what follows we will examine the poles of the Green's function G{p,y). The 
positions of these poles depend on the functions Ghom(p, 2/), and Dq as defined in 
Eqs. (21) and (22). The choice of a particular rule of treatment of the poles cor- 
responds to the choice of appropriate boundary conditions in the coordinate space. 
Note that the latter are dictated by physical constraints on the Green's function G 
rather than on the auxiliary function D. 

To get to the main point, we will try the simplest strategy of specifying the poles 
and check, aposteriori, whether this strategy is self-consistent. Let us put 

c = 

and define D in the Euclidean momentum space. Since in the Euclidean space the 
expression for D is well-defined and has no singularities, 

D{PE, q) = , D{pE, q)= I d^'y e'^^ D{pe, y) , (28) 

q'=j:ii''r: 

a 

one can perform analytic continuation from the Euclidean space to Minkowski. This 
is not the end of the story, however. It is the Green's function G that we are 
interested in, not the auxiliary function D. As will be explained below, the above 
procedure is consistent, for the following reason. The function G obtained in this 
way has a cut extending from zero to infinity. In addition, we find two complex 
conjugate poles on the second non-physical Riemann sheet of the complex p^ plane. 
Moreover, there are additional poles on subsequent unphysical sheets. 

Since the poles are not on the physical Riemann sheet, they do not correspond 
to any asymptotic states of the theory. A pole on the second Riemann sheet is 
a well-known signature of a resonance state [41]. Therefore, our toy scalar "grav- 
ity" is mediated by a massive resonance. The resonance-mediated gravity was first 
discussed in Refs. [42, 43, 44] in a different brane-world model. 

Before passing to consideration of particular cases it is worth reminding that the 
Green's functions in the N >2 DGP models need a UV regularization [2, 31, 32, 26]. 
This has been already mentioned. An appropriate UV regularization can be achieved 
either by introducing an explicit UV cutoff, or, alternatively, by keeping a non-zero 
brane width in a consistent manner (defined in Ref. [26]). For brevity we choose the 
former prescription by consistently taking the limit of zero brane width. However, we 
should stress that all our results hold equally well in the brane-width regularization 
method of [26]. 
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2.3 Six dimensions 



It is instructive to demonstrate how things work by considering separately the six- 
dimensional case. In six dimensions sensitivity to the UV cutoff is only logarithmic, 
and it is conceivable that the results obtained in the cut-off theory could be consis- 
tently matched to those of a more fundamental UV-completed theory-to-come. ^ 
It is not difficult to calculate 



(29) 



where is an ultra-violet cut-off. With this expression for Dq the function G{p'^, 0) 
develops a cut on the positive semi-axes of s due to the logarithmic behavior of 
Dq{s). This fact has a physical interpretation. Since the extra dimensions are non- 
compact in the model under consideration, the spectrum of the theory, as seen from 
the 4D standpoint, consists of an infinite gapless tower of the KK modes. This 
generates a cut in the Green's function for s ranging from zero to +cx3. 

In addition, there might exist isolated singular points in G(p^,0). These singu- 
larities (for s <^ A^) are determined by the equation 



G^\s, 0) = s 




1 -1 



0, 



(30) 



where is defined in Eq. (5). Let us introduce the notation 

s = So exp(i7) , (31) 
where Sq is a real positive number. Then, Eq. (30) has two solutions of the form 



So ~ 47r ml 



In ^ 



and 



7i - - 



TT 



log(A7m2) 



72 ~ 27r 



TT 



log(A2/m2) 



(32) 



(33) 



We conclude that there are two complex-conjugate poles on the nearby non-physical 
Riemann sheets. These poles cannot be identified with any physical states of the 
theory. They are, in fact, manifestations of a massive resonance state. All other 
complex poles appear on subsequent nonphysical Riemann sheets. 



^Thc D > 6 models of brane-induced gravity are power sensitive to UV physics. In general one 
expects all sorts of higher derivative operators in this case. 
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2.4 More than six dimensions 



Physics at D > 6 is similar to that of the six- dimensional world which was de- 
scribed in Sect. 2.3. There are minor technical differences between odd- and even- 
dimensional spaces, however, as we will discuss momentarily. 
In seven dimensions we find 

Do{s) = ^ |a - V^arctan [ ■ (34) 



As in the 6D case, there is a branch cut. The cut in this case is due to the depen- 
dence of the Green's function on No other singularities appear on the physical 
Riemann sheet. All poles are on non-physical Riemann sheets, as previously. 
In the eight-dimensional space the expression for Dq reads 



Again, we find a cut due to the logarithm, similar to that of the 6D case. All isolated 
singularities appear on non-physical Riemann sheets. 

The nine-dimensional formula runs parallel to that in seven dimensions, 

Do{s) = |y + s (a - V^arctan-^j | . (36) 



Finally, in ten dimensions 

1 



Do{s) 



128 7r3 




A^ + s Mn ^ + 1 



(37) 



The pole structure of G is identical to that of the eight-dimensional case. Since the 
pattern is now well established and clear-cut, there seems to be no need in dwelling 
on higher dimensions. 

Before turning to gravitons we would like to make comments concerning the 
UV cutoff A. The crossover distance Vc ~ depends on this scale: in 6D the 
dependence is logarithmic, while in D > 6 this dependence presents a power-law 
[2, 4]. Hence, the crossover scale in the N > 2 DGP models, unlike that in the 
N = 1 model, is sensitive to particular details of the UV completion of the theory. 
Since in the present work we adopt an affective low-energy field-theory strategy, we 
are bound to follow the least favorable scenario in which the cutoff and the bulk 
gravity scale coincide with each other and both are equal to Af* ~ 10~^ eV. If a 
particular UV completion were available, it could well happen that the UV cutoff 
and bulk gravity scale were different from the above estimate. In fact, in the string- 
theory-based construction of Ref. [4] the UV completion is such that the cutoff and 
bulk gravity scale are in the ballpark of TeV. 
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In conclusion of this section it is worth noting that the Green's function Dq in the 
N >3 case contains terms responsible for branch cuts. These terms are suppressed 
by powers of s/A, and, naively, could have been neglected. It is true, though, 
that the explicit form of these terms is UV-sensitive and cannot be estabhshed 
without the knowledge of UV physics. One should be aware of these terms since 
they reflect underlying physics — the presence of the infinite tower of the KK states. 
Fortunately, none of the results of the present work depend on these terms. 

3 The graviton propagator 

Now it is time to turn to gravitons with their specific tensorial structure. We will 
consider and analyze the equation of motion of the DGP-type model presented in 
Eq. (8), which we reproduce here again for convenience 

5W(l/)M^i (aG(f) + bGl^J)S^^S^^ + Ml^^G^^i = -T,^{x) 5^^5'^^5^^\y) . (38) 

Here G*^^-* and G^^^ denote the four-dimensional and D-dimensional Einstein tensors, 
respectively, while a and h are certain constants satisfying the constraint 

a + 6 = 1. 

For simplicity we choose a source term localized on the brane, namely, T^^(x)d^^\y). 
At the effective-theory level the ratio a/6 = a/(l — a)isa free parameter. The 
only guidelines we have for its determination are (i) phenomenological viability; 
(ii) intrinsic self-consistency of the effective theory which, by assumption, emerges 
as a low-energy limit of a self-consistent UV-completed underlying "prototheory." 
Specifying the prototheory would allow one to fix the ratio a/(l — a) in terms of 
fundamental parameters. 

Our task is to study the gravitational field produced by the source T^^{x)6^^\y). 
To this end we linearize Eq. (38). If qab = Vab + '^hAB, in the linearized in h 
approximation we find 

= dlhAB ~ dAdch% - dndch'^i 

+ Oa Ob /ig - VAB dl /ig + r)AB dc do /i^"" , (39) 

where d'jj = dod^. On the other hand, the four-dimensional Einstein tensor in the 
linearized approximation is 

- V>.^djh2 + r)^,dadph^r (40) 
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In what follows we will work in the harmonic gauge, 

d^hAB ^Idsh'^. (41) 

The advantage of this gauge is that in this gauge the expression for G^^ significantly 
simplifies, 

= dlhAB - IvABdlh^. (42) 

Additional conditions which are invoked to solve the {ab} and {a/i} components of 
the equations of motion are 

Ki^ = 0, hab = ^ Vab he . (43) 
Using the last equation it is not difficult to obtain the relation 

Nh^^^ {2-N)hl. (44) 

This relation obviously suggests that we should consider separately two cases: 

(i) N = 2; 

(ii) N>2. 

We will see, however, that the results in the N = 2 and N > 2 cases are somewhat 
similar. 



3.1 Brane-induced gravity in six dimensions (AT = 2) 

In two extra dimensions Eq. (44) implies 

h', = 0. (45) 

Therefore, the trace of the £)-dimensional graviton coincides with the trace of the 
extra-dimensional part, 

hi^K. (46) 

As a result, the four-dimensional components of the harmonic gauge condition (41) 
reduce to 

d^h^. = ld.K. (47) 

Let us now have a closer look at the {^i^} part of Eq. (8). Taking the trace of this 
equation and using Eqs. (42), (40), (45) and (47) we arrive at^ 

(36- l)<5W(y)Mg,a; K + 2 M,^+^ 9^ K = S^""^ iv) ■ W 
^As before, we put the transverse derivatives to be zero in the A — > hmit. 
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The obtained equation is very similar to the scalar- field equation (12). Therefore, 
we will follow the same route as in the scalar-field case, until we come to a subtle 
point, a would-be obstacle, which was non-existent in the scalar-field case. 
Let us Fourier-transform Eq. (48), 



{3b-l)6^''\y)M^,{-p')h:ip,y) 
+2Mr^ (V - A,) Kip, y) = f{p) <5W(y) . 
The general solution of the above equation is 

f{p) 



KM 



D{p,y) 



+ eg, 



horn : 



2m2 - (36 - l)p^Do{p) 
where the solution of the homogeneous equation takes the form 
e^hom = D{p,y)5[2ml - (36 - l)p'Do(p)) . 



(49) 

(50) 
(51) 

(52) 



To begin with, let us consider the case 36 > 1. Then the first term on the right-hand 
side of Eq. (51) has poles for complex values of p^, as can be readily seen from the 
expressions for Dq obtained in Sect. 2. For instance, in the 6D case this pole is 
determined by the equation 



2ni-: 



An 2ml 



(36 - l)A)(.s) (36 - 1 
This equation has at least two solutions of the form 



—s 



(53) 



47r 2ml 
36-1 



1 



and 



TT 



71 



log(A^/m- 



72 



27r + 



TT 



log(A2/rri 



(54) 



(55) 



The quantity h'^{p,y) is not a gauge invariant variable. Therefore, the presence of 
certain poles in the expression for h'^{p,y) depends on a gauge. However, explicit 
calculations (see below) show that the poles found above also enter the gauge in- 
variant physical amplitude. Therefore, we need to take these poles seriously and 
analyze their physical consequences. 
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3.2 6>l/3 

If 6 > 1/3 there are no poles on the physical Riemann sheet. Instead, poles appear 
on the nearest non-physical Riemann sheets. These poles cannot be identified with 
any physical states of the theory. They represent a signature of massive resonance 
states. All other complex poles appear on subsequent nonphysical Riemann sheets. 

Using a contour integral one can easily write down the spectral representation 
for the Green's function Q 

, I f>(t)dt , , 

TT ^0 t — 'p^ — le 

where the spectral function is defined as 

^ 2mglmZ}o(t) 

[(36- l)tReDo - 2m2]2 + [(36- l)tImDo)]^ ' 



and 



iV+2 



N-2 



= ^ y (2# = (2.rm/2) *^ ■ 

We see that p{t) satisfies the positivity requirement. Equation (56) guarantees that 
the Green's ftmction Q is causal. 

The next step is applying the expression for Q to calculate /i^i^. In fact, it is 
more convenient to calculate the tree-level amplitude 

A(p,y) = h,.{p,y)f^^(p), (59) 

where T'^^^[p) is a conserved energy momentum tensor, 

p^f'^"" = p^f'^"" = 0. 

Using Eqs. (38), (50) and (70) we obtain the following expression for the amplitude 
A{p,y): 



(26-l)p2Do 



ml 



(36 - 1)p^Dq - 2ml 



(60) 



Let us study the above expression in some detail. The first question to ask is about 
poles. It is quite clear that the p^-poles of A are of two types; their position is 
determined by: 

P^Dq{p) = ml 

or 

(36-l)p'L>o(p) = 2m^ 
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As was explained previously, all these poles appear on the second Riemann sheet, 
with the additional images on other non-physical sheets. None of these poles can 
be identified with asymptotic physical states. As was elucidated above, the occur- 
rence of the poles on the second and subsequent Riemann sheets corresponds to the 
massive-resonance nature of the effective 4D graviton. Our previous analysis can 
be repeated practically verbatim, with minor modifications, proving analyticity and 
causality of the amplitude A. 

Next, we observe that at large momenta, i.e., when Dq{p) ^ ml, the scalar 
part of the propagator has 4D behavior; the tensorial structure is not four-dimensional, 
however. The terms in the braces in Eq. (60), namely. 



correspond to the exchange of massive gravitons and scalar degrees of freedom. This 
would give rise to additional contributions in the light bending, and is excluded 
phenomenologically, unless the contribution due to extra polarizations is canceled 
by some other interactions (such as, e.g., an additional repulsive vector exchange). 
Note also that when b ^ a, i.e., 6 — > 1, one obtains the tensorial structure of 6D 
gravity, as expected from (7). 

On the other hand, at large distances, i.e., at Do{p) <^ ml, we get the following 
tensorial structure of the amplitude (60): 



This exactly corresponds to the exchange of a six-dimensional graviton, as was 
expected. 



This case is conceptually different from that of Sect. 3.2. As we will see momentarily, 
if 6 < 1/3 there are no problems in (i) maintaining 4D unitarity; and (ii) getting 
the appropriate 4D tensorial structure of gravity at sub-horizon distances. This is 
achieved at a price of abandoning 4D analyticity, in its standard form, which could 
presumably lead to the loss of causality at distances of the order of m~^ ~ 10^^ cm. 
The absence of causality at distances > 10^^ cm, was argued recently [27] to be an 
essential ingredient for solving the cosmological constant problem. 

Although all derivations and conclusions are quite similar for any ratio a/b as 
long as 2b < a, we will stick to the technically simplest example 6 = 0, a = 1. In 
the situation at hand, the homogeneous part (52) need not be trivial, i.e c need 
not vanish. The value of the constant c is determined once the rules for the pole at 
Dq{p) -\-2ml — are specified. In Ref. [40] the vanishing of c was postulated. This 




(61) 




(62) 



3.3 h < 1/3 
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choice leads to non-unitary Green's function. Therefore, wc abandon the condition 
c = in an attempt to make a more consistent choice that would guarantee 4D 
unitarity. We stress that we are after unitarity here, not unitarity plus causality. 

To begin with we pass to the Euchdean space in (i.e. p%) and introduce 
the following notation: 

^'^'(-y - 2r„i-M.)-i^ - '''' 

The function P^^^ is a Euclidean-space solution of Eq. (49), with the particular 
choice c = iir. (The choice c — —iir would lead to Eq. (63) with the replacement 
6-.-e). 

As the next step we will analyze the complex plane of p^. Since the function 
Dq{pe) is real, the function P'^^\p'e) must (and does) have an isolated singularity 
in the p\ plane which is similar to a conventional massive pole, except that it lies 
in the Euclidean domain. This singularity occurs at the point p\ — pi is defined by 
the condition 

pIDq^p^) = 2ml, P* real and positive. (64) 

This is the only isolated singularity in Eq. (63); it is located in the complex p^. plane 
on the real positive semiaxis. In addition to this pole singularity, the function (63) 
has a branch cut stretching from zero to — oo due to the imaginary part of Do(p|;) 
appearing at negative values of p%. As before, this branch cut is the reflection 
of an infinite gapless tower of the KK states. As a result, the following spectral 
representation obviously emerges for P^^\p%): 

^(F\, 1 1 r"^ l-m.P'^^Xu) du R 

p(^){pl + ze) ^ - 2 65 

n Jo u - p% pi - Pe - 

with the Euclidean pole term being "unconventional." The residue of the pole R is 
given (for any N) by 

J {2n)^ (g2 + p2)2 • 1^0^ 

Note that in the first term on the right-hand side of Eq. (65) the integration 
runs from zero to minus infinity; thus, the integrand never hits the would-be pole at 
u = p% > 0. Therefore, the ie prescription is in fact used only to specify the isolated 
pole at p%=pI. 

We proceed further and define a symmetric function 

n^^Hpl) = \ {p^^\pI - ie) + P^^^ipl + ^6)} . (67) 



18 



It is just this symmetric function on which we will focus in the remainder of the 
section. Let us return to Minkowski space. This is done by substituting 

p% — > exp(— i7r)p^ , u — > exp{—in)t 

in Eq. (65). Furthermore, observing that ImP — Imll, we obtain the following 
representation for the Minkowskian H: 

1 lmU(t)dt 

n p = - / + Ho p , 68 

TT Jo t — — le 



where 



no(p) = \ — - + ^-4— -1 ■ (69) 

2 \pi - le pi + lej 

It is necessary to emphasize that e and e are two distinct regularizing parameters, 
e 7^ e. The parameter e is used to regularize the pole at p^ — —pi, while e sets 
the rules for the branch cut. The most important property of 11 is that the pole 
at p"^ = —pi has no imaginary part, by construction. Hence, there is no physical 
particle that corresponds to this pole. In the conventional local field theory the only 
possible additions with no imaginary part are polynomials. Here we encounter a 
new structure which will be discussed in more detail at the end of this section. 

Our goal is to show that a 4D-unitarity-compliant spectral representation holds 
for the Green's function on the brane, at least in the domain where the laws of 4D 
physics are applicable. To this end we turn to the function Qip,y), defined as 

g = D{p,y)U{p'). (70) 

with the purpose of studying its properties. It is convenient to pass to the momentum 
space with respect to extra coordinates too. Then, the propagator (70) takes the 
form 

G{p,q) = , ■- ■ (71) 

— p^ — le 

With these definitions in hand, we can write down the 4D dispersion relation. We 
start from the Kallen-Lehman representation for the propagator (71). As we will 
check below, this representation takes the form 

gfa,)^ir + (72) 

TT JO t — p^ — le q'^ — p'^ — le 
The imaginary part of Q is defined as follows 

Im^(i, q) = 7r(5(g2 - t) Reli{t) + Imn(i) , (73) 
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where V stands for the principle value of a singular function, 



g2 - t 




(74) 



The fact that Eq. (72) holds can be checked by substituting (74) and (73) into (72) 
and exploiting the relation 

1 r Imn(t) ^ 1 Ren(g^) - no(g^) + Yi,{p^) - Ii{p^) 

~ / 1 2 ~'^~2 7"* ^ 2 2 ~ ■ y'^' 

TT Jo t — p"^ — te q"^ — t q^^ — p^ — le 

This turns Eq. (72) into identity. 

Finally we approach the main point of this section - the dispersion relation for 
S{p,y = 0), the Green's function on the brane. As such, it must have a spectral 
representation with the positive spectral density, as we have already seen from the 
KK-based analysis. The positivity is in one-to-one correspondence with the 4D 
unitarity. 

The dispersion relation can be obtained by integrating (72) with respect to q, 

Gip, y = 0) = - r ^ P^Y* ._ + Mp^) ReD,{-pl) . (76) 
n Jo t — p'^ — le 

According to Eq. (72), the spectral density p is defined as 



Pit) = I (77) 



The first term on the right-hand side in Eq. (76) is conventional while the second is 
not, and we hasten to discuss it. This term has no imaginary part, by construction. 
Hence, it does not contribute to the unitarity cuts in diagrams. Therefore, this term 
does not affect the spectral properties. 

As was mentioned, in conventional 4D theories only a finite-order polynomial in 
that has no imaginary part can be added to or subtracted from the dispersion 
relation. This is because, normally one deals with Lagrangians which contain only 
a finite number of derivatives, i.e., a finite number of terms with positive powers 
of p^ in the momentum space. In the problem under consideration this is not the 
case, however. In fact, no local 4D Lagrangian exists in our model at all, and yet 
we are studying the spectral properties in terms of the intrinsically 4D variable, 
p^. The theory (1) is inherently higher-dimensional because of the infinite volume 
of the extra space. One can try to "squeeze" it in four dimensions at a price of 
having an infinite number of 4D fields. For such a theory there is no guarantee that 
analyticity of the Green's functions in terms of the 4D variable p^ will hold because 
the effective 4D Lagrangian obtained by "integrating out" the infinite gapless KK 
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tower will necessarily contain [5] non-local terms of the type d~'^. (Note that a 
similar prescription for the poles in a pure 4D local theory [45] is hard to reconcile 
with the path integral formulation [46]. In our case this is not a concern since the 
theory is not local in four-dimensions in the first place.) 

Therefore, it is only natural that 4D unitarity can be maintained h\ii 4D ana- 
lyticity cannot. Non-analyticity leads to violation of causality, generally speaking. 
That is to say, the Green's function (76) is acausal. Therefore, we have an apparent 
violation of causahty in the 4D slice of the entire (4 + N) dimensional theory which, 
by itself, is causal. The apparent acausal effects can manifest themselves only at the 
scale of the order of ~ 10^^ cm. In fact, as was noted in [27], this is a welcome 
feature for a possible solution of the cosmological constant problem. 

Let us now return to the first term on the right-hand side of Eqs. (76). Using 
Eqs. (68) and (73) we can calculate the spectral function which comes out as follows: 

n(f\ = 2mglmDo(t) 

' {tReDo + 2m2)2 + (tImDo)' ' 

where 

We see that p(t) satisfies the positivity requirement. ^ 

Next we observe that at large momenta, i.e., at Dq{p) ^ ml, the propagator 
we got has the desired 4D behavior. For the scalar part of the propagator this is 
expected from the studies of Sect. 2. However, with regards to the tcnsorial structure 
this circumstance is less trivial. If Dq{p) ^ ml the terms in the braces in Eq. (60), 

f^uf'^" - f'l , (80) 

correspond to the exchange of two physical graviton polarizations. Therefore, for 
the observable distances the tensorial structure of the massless 4D graviton (80) is 
recovered. 

On the other hand, for large (super- horizon) distances, p^ Dq{p) <C m^, we get a 
different tensorial structure of the same amphtude, 

If^.f^'' -]ff' \ . (81) 



4 

This exactly corresponds to the exchange of the six- dimensional graviton 



^For iV > 5 the integral in Eq. (76) diverges. However, since our model has a manifest UV 
cutoff A, the above integral must be cut off at A. Alternatively, one could use a dispersion relation 
with subtractions. 
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3.4 D > 6 

Corresponding calculations and results are quite similar to the D = 6 case, with 
minor technical distinctions which we summarize below. For N ^ 2 

Kb^J^VabK- (82) 



Therefore, we get 



d'K^- J^d^K- (83) 



Then, the trace of the {liv} equation takes the form 

kN 5^""^ (y) Ml, dl K + {N + 2) d\ K 

= NTI^S^'^^y). (84) 

where kN = 2 — iV(2 — 36). The above equation which can be used to find the solution 
we are after. We proceed parallel to the six-dimensional case. Let us introduce the 
notation 

kip,y) = ^^^^^M, (85) 

where 

^ D{P, y) , r (Qa\ 

-KnP Do{p) + [N + 2) 

The solution of the homogeneous equation takes the form 

^^Nhom = D{p, y) S {-k^p' Do{p) + «^ (iV + 2)) . (87) 

Here 

As in the 6D case, we conclude that that there exists a solution to the equation 

-kNP^Doip) + u^{N + 2) = 

with a complex value of p^. These poles occurs on the nonphysical sheets as long as 
kjsf > 0, so the Green's function admits the spectral representation. 

Using the expressions above one readily calculates the tree-level amplitude A, 

1 Dip,y) 



^^^'^^ Ml, p^ Doip) - 

I kNP^Do-{2 + N)u^ ^^^^ 
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3.5 6> (2iV-2)/3iV 



In this case there are no poles on the physical Riemann sheet. Hence, all the poles 
are of the resonance type. The tensorial structure at large distances is that of the 
D-dimensional theory 

[t^.T'^'^ - ^fry (89) 

However, the tensorial structure at short distances ^ differs from that of 4D 
massless gravity. Hence, some additional interactions, e.g., repulsion due to a vector 
field, is needed to make this theory consistent with data. 



3.6 6 < (2iV- 2)/3iV 

The consideration below is very similar to the 6D case. In perfect parallel with the 
6D case we consider for simplicity only the 6 = case and define the function 

P^n\pI) = uN{N + 2)/{2N-2) -plDo{pE)-ie' 



which has a spectral representation: 

u - Pe Pi - Pe - 

The residue R is determined by Eq. (66) while pi is now a solution to the equation 



p'/u = - r 5=^^^ + , 1 . . (91) 

TT Jo U - pi, Pt - Ptt. - le 



PlDoiP.) - "^^^^ , pI>0. (92) 
As in the 6D situation, we use the expression (91) to define a symmetric function 

<\pI) ^ \ {p^n\pI - ie) + HpI + ie)}. (93) 

The latter, being continued to the Minkowski space, admits the following spectral 
representation: 

n„(rt = i r-M^_ + nnp). (94) 

TT Jo t — p^ — le 

where 

nr(,)4(^;^^ + ^;^^). (95) 
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As previously, e and e are two distinct regularizing parameters, e ^ e. 
For the Green's function of interest 



D{p,y) n^(p2) 



(96) 



2N -2 



we repeat the analysis of Sects. 3.1, 3.2 and 3.3 to confirm with certainty that the 
function QN{p,y = 0) does admit the spectral representation (76), with a positive 
spectral function, similar to the 6D case, see Eq. (78). 

The expression in Eq. (88) interpolates between the four-dimensional and D- 
dimensional patterns. This was already established for the scalar part of the ampli- 
tude in Sect. 2. Let us examine the tensorial part. For -Do(p) ^ we get 



This corresponds to two hehcities of the 4D massless graviton. In the opposite hmit, 
Do{p) <^ , we recover the tensorial structure corresponding to the (4 -|- A^)- 
dimensional massless graviton. 

4 Perturbation theory in massive gravity: 
hard mass vs. soft 

We start from a brief review of the well-known phenomenon — the breakdown of 
perturbation theory for the graviton with the hard mass [34] , occurring at the scale 
lower that the UV cutoff of the theory [35, 36, 37]. We then elucidate as to how this 
problem is avoided in the models (1, 7). 
The 4D action of a massive graviton is 



Sm-^l d^x^gRig) + I d'x {hi - {hff) , (98) 



where nig stands for the graviton mass and h^i, = {g^i, —r]n^)/2. The mass term has 
the Pauli-Fierz form [34] . This is the only Lorentz invariant form of the mass term 
which in the quadratic order in h^i, does not give rise to ghosts [47] . Higher powers 
in h could be arbitrarily added to the mass term since there is no principle, such as 
reparametrization invariance, which could fix these terms. Hence, for definitcness, 
we assume that the indices in the mass term are raised and lowered by rj^^. Had we 
used Qfj^i, instead, the difference would appear only in cubic and higher orders in h, 
which are not fixed anyway. 

In order to reveal the origin of the problem let us have a closer look at the free 
graviton propagators in the massless and massive theory. For the massless graviton 




(97) 
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we find 



Dlu;al3iP) = Vf^aVup + ^ Vt^pVua ' \ Vf^uVaP^ _ _ , (99) 

where 

f\ixv = Vi^i^ - ■ (100) 

The momentum dependent parts of the tensor structure were chosen in a particular 
gauge convenient for our discussion. On the other hand, there is no gauge freedom for 
the massive gravity presented by the action (98); hence the corresponding propagator 
is unambiguously determined, 



where 



V,. ^V,.-^. (102) 

mi 



We draw the reader's attention to the I /nig, singularities of the above propa- 

gator. The fact of their occurrence will be important in what follows. 

It is the difference in the numerical coefficients in front of the rj^i^rjap structure 
in the massless vs. massive propagators (1/2 versus 1/3) that leads to the famous 
perturbative discontinuity [48, 49, 50]. No matter how small the graviton mass is, 
the predictions are substantially different in the two cases. The structure (101) gives 
rise to contradictions with observations. 

However, as was first pointed out in Ref. [35], this discontinuity could be an 
artifact of relying on the tree-level perturbation theory which, in fact, badly breaks 
down at a higher nonlinear level [35, 36]. One should note that the discontinuity 
does no appear on curved backgrounds [51, 52] — another indication of the spurious 
nature of the "mass discontinuity phenomenon." 

To see the failure of the perturbative expansion in the Newton constant Gn 
one could examine the Schwarzschild solution of the model (98), as was done in 
Ref. [35] (sec also [53, 54, 55]). However, probably the easiest way to understand 
the perturbation theory breakdown is through examination of the tree-level trilinear 
graviton vertex diagram. At the nonlinear level we have two extra propagators which 
could provide a singularity in nig up to l/rrig. 

Two leading terms, l/rrig and I /nig, do not contribute [36], so that the worst 
singularity is l/m^- This is enough to lead to the perturbation theory breakdown. 
For a Schwarzschild source of mass M the breakdown happens [35, 36] at the scale 

A^~m,(Mm,/M2i)-V5. 
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The result can also be understood in terms of interactions of longitudinal polariza- 
tions of the massive graviton which become strong [37]. For the gravitational sector 
per se, the corresponding scale reduces to [37] 



;m,/MpO-^/^ 



If one uses the freedom associated with possible addition of higher nonlinear terms, 
one can make [37] the breaking scale as large as 

m,/K/MpO^/^ 

(Note that at the classical level the strong-coupling problem of the PF gravity can be 
evaded by summing up tree- level nonlinear diagrams [35, 36]. To determine whether 
the problem is present at the quantum level, one must perform perturbations on a 
stable background; however, the Minkowski-space background is not stable for the 
PF gravity, with the instability setting in almost instantaneously [39]. For recent 
discussions of massive gravity see Refs. [56], [57].). 

Summarizing, in the diagrammatic language the reason for the precocious break- 
down of perturbation theory can be traced back to the infrared terms in the propa- 
gator (101) which scale as 

(103) 

These terms do not manifest themselves at the linear level; however, they do con- 
tribute to nonlinear vertices creating problems in the perturbative treatment of 

massive gravity already in a classical theory. 

We will see momentarily that similar problems arc totally absent in the propa- 
gator of the model (1). For illustrational purposes it is sufficient to treat the N = 2 
case. All necessary calculations were carried out in Sect. 3. Therefore, here we just 
assemble relevant answers. 

For AT = 2 and b > 1/3 we find 

P,P:^D{p.y) ^^^^^ 



2ml - (36- l)p2Do(p) + ie 



In the limit mc — > the above expression, as opposed to Eq. (103), is regular. Similar 
calculations can be done in the N > 2 case. The results is proportional to 

P,P.D{p,y) ^^^^^ 



{2 + N)u^ - kNp'^Doip) + ie' 



which is also regular in the ttIc limit where it approaches the 4D expression. 
Therefore, we conclude that there is no reason to expect any breaking of perturbation 
theory in the model (1) below the scale of its UV cutoff. 
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If 6 < 1/3 and N = 2 we find, by the same token, 



Again, in the hmit mc — > the above expression, in contradistinction with Eq. (103), 
is regular. Moreover, in this hmit (and at y = 0) it approaches the 4D expression, 
in a particular gauge. Analogous calculations can be readily done in the N > 2 and 
b < {2N - 2)/3N case. The results is 



1 

+ (e ^ -e) 



{2 + N)u^ -kNp" Do{p)+ie 



(107) 



This expression is also regular in the mc — > limit where it arrives at the correct 4D 
limit. We conclude therefore, that in the general case there is no reason to expect 
any breaking of perturbation theory in the model (7) below the scale of its UV 
cutoff. Note that the expressions (106) and (107) are singular for small Euclidean 
momenta ~ ~^c- construction this singularity has no imaginary part and 
there is no physical state associated with it. One might expect that this singularities 
will be removed after the loop corrections are taken into account in a full quantum 
theory. These considerations are beyond the scope of the present work. 

An analogy with the Higgs mechanism for non-Abelian gauge fields is in order 
here. For massive non-Abelian gauge fields nonlinear amplitudes violate the unitar- 
ity bounds at the scale set by the gauge field mass. This disaster is cured through 
the introduction of the Higgs field. Likewise, nonlinear amplitudes of the 4D massive 
gravity (98) blow up at the scale Am- The unwanted explosion is canceled at the 
expense of introducing an infinite number of the Kaluza-Klein fields in (1). 



5 Discussion and conclusions 

In the present work we studied the model (7) of the brane-induced gravity in codi- 
mensions two and higher. This model has stringent and testable predictions. Gravity 
is modified at short distances, of the order of ~ 0.1 mm or so, and simultaneously, 
at ultra-large distances, of the order of ~ 10^^ mm, give or take an order of mag- 
nitude. The short-distance modification can be tested in table-top gravitational 
experiments [58, 59]. Modification of gravity at a millimeter scale and its relation 
to the cosmological constant problem was first discussed in Ref . [60] . 

The modification of gravity at a millimeter scale in the present model (1) is a 
consequence of the large-distance modification and vise versa. These are two faces 
of one and the same phenomenon. However, we should point out that technically 
and conceptually the approaches to the cosmological constant problem discussed in 
Ref. [60] on the one hand, and Refs. [2, 5] on the other, are rather different — the 
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former relies on the short-distance (UV) modification of gravity, while the latter is 
entirely based on the large- distance (IR) modification of gravity. 

The large distance modification of gravity can manifest itself in cosmological so- 
lutions. The case 6 = and a — 1 seems to be most interesting for these purposes. 
As we argued in the present work, it leads to apparent violations of 4D causality 
which could manifest themselves at the scales of the order of today's Hubble scale. 
Manifestations of acausality might be tested in cosmological observations. In par- 
ticular, such an acausal theory might be the reason behind the smallness of the 
observable space-time curvature [27] . 

It is instructive to point out how the 6 = 0, a = 1 model evades a well-known 
no-go theorem for massive gravity [49]. Let us first bricfiy recall the theorem. A 
4D massive graviton has extra polarizations one of which couples to sources in the 
leading order in a weak field. The additional attraction due to this polarization 
is observationally unacceptable and has to be canceled. This can be achieved by 
introducing a ghost that gives rise to compensating repulsive force [49]. Hence, either 
one ends up with a theory that has a ghost or with a theory that has no ghosts but 
is phenomenologically unacceptable. This is the essence of the no-go theorem [49]. 
The theorem can easily be generalized for a theory with an infinite number of states 
[61]. In the latter case, in order to obtain a phenomenologically acceptable theory 
of a massive graviton at observable distances one should give up positivity of the 
spectral function in the dispersion relation for the corresponding Green's function. 
This would violate unitarity of the model. However, the above argument assumes 
4D analyticity the consequence of which is the existence of a spectral representation 
for the Green's function. In our case 4D analyticity is violated, and so are the 
conditions of the no-go theorem. 

We would like to point out a certain common feature with (2 + 1) topologically 
massive gauge/gravity theory [62] where the large distance interactions are also 
power-like. 

Finally, we would like to emphasize that the models (1) and (7) give rise to a 
gravitational analog of the Higgs mechanism in the following sense — the effective 
graviton-mediating interaction is massive, and, nevertheless, the growth of the non- 
linear amplitudes is softened at the expense of having an infinite number of fields. 
This phenomenon can be rather transparently understood from the standpoint of 
the KK modes. The manifest reparametrization invariance is a convenient book- 
keeping tool in this case for determining whether or not the amplitudes blow up. 
The reparametrization invariance at each KK level is maintained at the same KK 
level only in the linearized approximation. Nonlinear effects mix distinct KK levels 
under the reparametrization transformations [63], [64]. Hence, if the KK tower is 
truncated at some finite level, the breakdown of perturbation theory in nonlinear 
diagrams is inevitable. However, if the infinite totality of the KK modes are kept, 
as in Eq. (1), the softening of the amplitudes should be expected. The present work 
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fully confirms this expectation. In light of this finding, it would be interesting to 
discuss the strong coupling issue in nonlinear interactions of the 5D DGP model 
studied in Refs. [36, 65, 66]. In particular, the immediate task is to understand 
whether this is a problem of peculiar perturbation theory, as advocated in Ref . [36] , 
or the problem inherent to the model itself [65, 66]. Already from the equations of 
the present work it is clear that the = 1 DGP model with 6 = 1 in (7) has no 
strong coupling problem. The other possibilities are currently under investigation; 
the answers will be reported elsewhere. 
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